Signatures of superfluidity for Feshbach-resonant Fermi gases 
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We consider atomic Fermi gases where Feshbach resonances can be used to continuously tune the 
system from weak to strong interaction regime, allowing to scan the whole BCS-BEC crossover. We 
show how a probing field transferring atoms out of the superfiuid can be used to detect the onset 
of the superfiuid transition in the high- Tc and BCS regimes. The number of transferred atoms, as 
a function of the energy given by the probing field, peaks at the gap energy. The shape of the 
peak is asymmetric due to the single particle excitation gap. Since the excitation gap includes 
also a pseudogap contribution, the asymmetry alone is not a signature of superfluidity. Incoherent 
nature of the non-condensed pairs leads to broadening of the peak. The pseudogap and therefore 
the broadening decay below the critical temperature, causing a drastic increase in the asymmetry. 
This provides a signature of the transition. 



Novel type of superfiuid, a Bose- Einstein condensate of 
molecules consistingof fermionic alkali atoms, has been 
created recently 0, 0]- Molecules are formed by a Fes- 
hbach resonance in the range of parameters where the 
interaction between the atoms is repulsive. By tuning 
magnetic fields the resonance can be crossed and attrac- 
tive interactions realized. Crossing of the resonance can 
be done adiabaticallv, in temperatures lower than the 
critical temperature [2j , and fcrmion pairing investigated 
01 . This opens up unprecedented prospects for studying 
the BCS-BEC crossover from a superfiuid of momentum- 
space pairs to a BEC of real-space dimers. The exper- 
iments 0, demonstrate the existence of a superfiuid 
in the BEC limit of the crossover. In this letter we pro- 
pose a method to observe an explicit signature of the 
superfiuid transition also in the intermediate (high-Tc) 
and BCS regimes. 

Fermionic alkali atoms in the vicinity of a Feshbach res- 
onance are predicted to form superfluids with three com- 
ponents: 1) fermionic Cooper-pairs, 2) molecular bosons 
formed out of the fermions, 3) non-condensed fermionic 
pairs corresponding to a pseudogap. Non-condensed 
pairs and molecules exist also above the critical temper- 
ature. The system can be described by the resonance 
superfluidity theory |^. In the BEC limit, the super- 
fluid transition is evident like for an atomic BEC, for 
instance via the bimodal density distribution P, |3| ■ In 
the high-Tc and BCS regimes such a strong effect on the 
density is not expected. Moreover, strong pseudogap or 
mean field may be almost indistinguishable from the ef- 
fects of superfluidity. Several signatures of superfluidity 
for fermionic atoms have been proposed in 7, 8], however, 
these works assume simple BCS scenario and do not take 
into account the pseudogap and the bosonic molecules. 
Here we discuss a method that can give a clear signature 
of the resonance superfluidity transition. It is experi- 
mentally feasible since a similar technique has already 
been used to study mean-field effects The method 
is based on transferring atoms from the superfiuid state 
into a normal state using electromagnetic fields. There 



is an analogy to imposing current over a superconductor- 
normal metal interface. Instead of a localized tunneling 
barrier, the transfer is caused by a homogeneous field, 
enforcing exact momentum conservation. This leads to 
a resonance-type of response as a function of the applied 
energy, in contrast to superconductors where the current 
grows steadily as a function of voltage once the gap en- 
ergy is exceeded . 

Interacting gas of two fermionic components can be 
formed by trapping atoms in two different hyperfine 
states, let us label the states by \g), \g'). The strength of 
the interaction between atoms in \g) and \g') is tuned by 
magnetic (or optical) fields using a Feshbach resonance. 
The atoms provide also other hyperfine states (label |e)) 
for which the interatomic scattering lengths can be very 
different from that between \g) and \g'). This is utilized 
in 9] to observe mean field effects: atoms are resonantly 
transferred, say, between \g) and |e}, and the difference in 
the mean fields experienced by \g) and |e) causes an en- 
ergy shift in the position of the resonance. Similarly, for 
a superfiuid formed of atoms in the \g) and \g') states, the 
BCS-gap in the single particle excitation spectrum is pre- 
dicted to cause a shift in the position of the resonance j3| . 
The main qualitative difference to the mean-field shift is 
the asymmetric shape of the peak in the superfiuid case. 
Here we show that a pseudogap can lead to a shifted, 
asymmetric peak as well. Therefore asymmetry alone is 
not a signature of superfluidity. The transition, however, 
is signaled by an abrupt change in the amount of the 
asymmetry. 

The system is described by the Hamiltonian 
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The operators correspond to fermionic atoms in var- 
ious hyperfine states {a — {e,^,^'}) and bg to spinless 
bosonic molecules formed out of atoms in states \g) and 
\g'). The fermion and boson energies are Sk = fc^/2m 
and = q'^/2M, respectively, with m the fermion 
mass, M = 2rn and h = 1. A momentum cutoff is 
used: U{k,k') = Uipkifik', g{k) = g(pk, where (fl = 
exp —{k/Kc)"^. The effect of the Feshbach resonance can 
be incorporated to a point-like interaction between the 
atoms, as in H, by using renormalized parameters U, 
g and v |^. The corresponding unrenormalized param- 
eters are: Uq = 47ras/m is the fermion- fermion inter- 
action, with as the scattering length, go is the boson- 
fermion coupling, and z/q the magnetic field detuning 
from the Feshbach resonance. Positive vq corresponds 
to attractive interactions between the atoms, negative 
Vq to repulsive. At the resonance, Vq = 0, the sys- 
tem is in the high- Tc region. The limits I'o — > oo and 
i>o — > — oo are the extreme BCS and BEC limits, respec- 
tively. The parameters used throughout the paper are 
go = —b.OEp and Uo = 0.05 Ep, corresponding to the 
density of atoms n = 10^^ cm"'^ and background scatter- 
ing length as — 176 cq. The initial numbers of atoms in 
states \g) and \g') are equal, therefore the correspond- 
ing chemical potentials /ig = ^g> = /i. The state |e) is 
assumed initially empty, that is, fie = 0. 

The observable of interest is the current of atoms from 
the state \g) to state |e), / =< >. The effect of the 
electromagnetic field transferring atoms out of the super- 
fluid, i.e. creating the superfluid-normal state interface, 
is described by Ht- To match the energy difference be- 
tween the considered hyperfine states, the field can be 
either an rf-field or laser field creating a Raman transi- 
tion. In describing the interaction of the field with the 
molecule (formed of fermions in states \g) and \g')), we 
have considered only the interaction of the field with the 
two- level system \g) — |e). This approximation is made 
having in mind large, extremely weakly bound molecules 
P,|3j|- We refer to the detuning S = {uj^ —ojg) —lo field by 
probe field detuning, not to confuse it with the Feshbach 
detuning vo of the magnetic field. 

The particle current I =< >^ i < [H,Nc] > is 
calculated treating Ht as a perturbation and neglecting 
terms of higher than second order. The solutions for the 
equilibrium parameters, chemical potential and the ex- 
citation gap, are needed for calculating the current. We 
apply the t-matrix approach to describe the pairing inter- 
action, adopting the procedure used in Ref. jllj for solv- 
ing the relevant equations: gap, pseudogap and number 
equations. The chemical potential ^, pseudogap Apg and 
the superfluid order parameter Age are obtained from the 



self-consistent solutions. 

Physically, one may interpret the pseudogap as fermion 
pairs which have a finite net momentum and are formed 
above the critical temperature (thus the name preformed 
pairs). At the critical temperature, these pairs condense 
to the zero net-momentum state. The origin of the pseu- 
dogap in the many-body theory are fiuctuations of the 
finite momentum order parameter (c^j.c^^^) around the 

zero-momentum mean (c^j.c^ ). The zero-momentum 
mean is the usual BCS order parameter, and has a non- 
zero value below the transition temperature. For finite 
momenta q, the mean {c^_i.cf._^^) is always zero. However, 
it has a non-zero variance, corresponding to fluctuations, 
also above Tc. This leads to the pseudogap. The i-matrix 
for the system can be written as a sum of a superfluid 
term A^^, containing a pole at Q — 0, and a pseudo- 
gap contribution that includes the rest of the t-matrix, 
Apg = ~X]q5^o*(*3) E3- The total excitation gap is 
given by A^ = A^^ -I- A^g . Pseudogap has a considerable 
effect in the high-Tc regime but vanishes in the extreme 
BCS limit. 

The pseudogap and the order parameter are qualita- 
tively different and in principle can be distinguished from 
each other through life-time effects. The incoherent na- 
ture of the propagating, non-condensed pairs affects scat- 
tering processes and thus incorporates a flnite broaden- 
ing in the self-energy Spg. This is taken into account by 
replacing the self-energy by the model self-energy ^l^ 

iLU + efc iLU + Ck + «7 

where 7 is the finite broadening. Using the model self- 
energy, we have plotted the fermionic density of states 
N{lu) = -2Y^f.lmG{uj + iO,k) in the inset of Fig. [H 
for several temperatures. We have chosen the Feshbach 
detuning vo = 2.0 iJp (attractive interactions) and the 
broadening 7 — A(Tc) « 0.6 i^p- The plot shows the 
widening of the peak at — —A as the temperature 
increases, caused mainly by the emergence of the pseu- 
dogap. At lower temperatures, a peak appears also at 
Lo = -|-A, and in the limit T ^ the density of states 
approaches the BCS limit of two peaks. 

The current I =< Nc >^ i < [H, Nc] > can now be 
calculated. Denoting A{t) = J^ki ^'^kic^ and B{t) — 

J2qki Dqkiblcf c|, where M^;, Dgki are the couphng con- 
stants, the correlations of the form < AA >, < BB >, 
< BA > and < B^A > (and h.c.) are neglected. The 
first three contain correlations of the form (c'^c'^) which 
are identically zero for a normal state. The last one cor- 
responds to a Josephson current between the molecular 
and fermion superfiuids, mediated by the |e) state [Hj . 

Writing the total current as a sum of 'fermionic' and 
'bosonic' currents / = /p + /b, one obtains using Mat- 
subara Green's functions techniques 
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where xf^ are the poles of the dressed Green's function 
Gg(n, z) = (z — — E(n, z))"^ and the coefficients ao 
and ai are obtained from the Unearisation of the boson 
propagator D{Q)~^ w ao{iuj — q^ai). The self-energy S 
receives an imaginary part from the broadening 7, giving 
also a nonvanishing imaginary part to the poles For 
7 = 0, the fermionic current If simplifies to the standard 
Fermi Golden rule form. We have neglected the backward 
process in the boson term /b as the |e)-states are initially 
unoccupied. The detunings are 5f = Sb = —fJ,+S—AE^{, 
where Ai^mf — — E^_^f is the difference in the mean 
field energies of the states \g) and |e). 

The current / as a function of the probe field detun- 
ing S is plotted in Fig. ^ The chosen Feshbach detun- 
ing i/Q = 2.0 Ep corresponds to the BCS side of the 
resonance (in the high-Tc regime, Tc w 0.27 i^p)- Due 
to small number of molecules on the BCS side, we use 
I — If + Ib ^ If- For comparison, we plot also the 
current neglecting the pseudogap (Apg = 0) in the lower 
figure. The main observations from the results are: 1) 
Effect of mean field interactions is simply a shift of the 
resonance peak by Ai^mf . This is shown by the symmet- 
ric Lorentzian at the center (solid line). It corresponds 
to a high temperature (~ Tp) with no preformed pairs. 
(Comparison: for Ap^ = 0, the simple mean field be- 
haviour persists until the critical temperature, the solid 
line Lorentzian is for T = Tc.) 2) Effect of pairing and 
excitation gap appears as asymmetric form and splitting 
of the peak. Pair formation occurs already above the crit- 
ical temperature, the dashed line for T = represents 
also T > Tc. Splitting of the peak is due to opening 
of an excitation gap. The side-peak on the left corre- 
sponds to breaking a pair: the gap energy A has to be 
provided, causing a shift from the center. The side-peak 
on the right originates from finite temperature quasipar- 
ticles, which are excitations with the minimum energy A. 
Energy is gained when transferring them from \g) to |e), 
therefore this side-peak appears on the opposite side of 
the center position. The quasiparticle side-peak is weaker 
and vanishes with decreasing temperature, leading to an 
asymmetric form of the total curve. An important fea- 
ture is that the curve is broadened by the finite lifetime 
of the pseudogap. (Comparison: for Apg = 0, pair forma- 
tion occurs only below Tc. Moreover, the pseudogap life- 
time broadening is absent which leads to sharper features. 
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Figure 1: Tlie particle current as a function of tlie probe 
field detuning 5 for different temperatures. The decay rate of 
the non-condensed pairs is chosen 7 = A(Tt;) ~ 0.6 -Ep. For 
comparison, the lower plot shows the current neglecting the 
pseudogap. The inset shows the density of states of fermions 
at resonance detuning uq = 2.0 -Ep. There are two peaks 
located at w = ±A, but increasing temperature broadens the 
peaks. The three plots are for temperatures T = Tc (solid), 
0.75 Tc (dashed) and 0.10 Tc (dotted). 



e.g. the two side-peaks are more visible. For Apg — 0, 
the existence of the excitation gap equals superfluidity.) 
3) Superfluidity is signaled by increase in the asymmetry 
of the peak (c.f. the dotted curve at T = 0.5 Tc) as will 
be discussed in the following. 

We quantify the asymmetry of the peak A = by 

''bright 

integrating the current / to the left wicft and to the right 
Wright from the peak value. The asymmetry has been plot- 
ted as a function of temperature in Fig. [3 signaling the 
phase transition by an abrupt change in the asymmetry 
as a function of temperature. The data for temperatures 
above Tc have been obtained by extrapolating the total 
excitation gap A(T) and solving the chemical potential 
fi from the number equation. Since the magnitude of 
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Figure 2: The asymmetry A of the current 7 as a function 
of temperature T/Tc for the Feshbach detuning vo — 2.0 Ef 
(attractive interaction). Above critical temperature (in the 
pseudogap regime) the asymmetry is essentially constant but 
rapidly increases in the superfluid regime. The lower solid 
line is for lifetime broadening 7 ~ A(r = Tc) and the upper 
line for 7 ~ 0.2 ■ A(r — Tc). For larger 7, the pseudogap 
pairs decay faster and are easier to distinguish from Cooper 
pairs, the lower line therefore shows sharper features at Tc- 
The dashed line neglects the pseudogap, showing the BCS- 
like behaviour. In the limit T —> 0, the pseudogap vanishes 
and all three curves approach the same limit. Averaging with 
relative densities n = 1, 0.8 and 0.6 is included. 



the lifetime broadening 7 is important for the effective- 
ness of the method, we shovif the asymmetry for two 7's 
of different orders of magnitude; 7 ^ A{T — Tc) and 
7 ~ 0.2A(r = Tc) in Fig. |2| Time-energy uncertainty 
relation suggests 7 — t^^^^y ~ A which is also used in 
[l5| . Note that the change in the asymmetry at T = Tc 
is visible also for the smaller 7. 

The main experimental challenge related to the 
method is likely to be mean field effects. Figure 
presents results for a homogeneous system where the 
mean field causes merely a shift of the resonance position. 
In real experiments, the gas is trapped in a harmonic po- 
tential and consequently the mean field shift varies de- 
pending on the local density in the trap. This leads to 
broadening of the shifted resonance peak Q- In order 
to estimate trapping effects, in the sense of local den- 
sity approximation, we have calculated the current and 
the asymmetry for relative densities n = 1.0, 0.8 and 0.6 
and mean fields AE^{{n) — 0.5 Ep. Fig. [21 presents av- 
eraged results: the change in the asymmetry is neverthe- 
less clear. The choice of the state |e) in such a way that 
mean field effects are under control, possibly combined 
with schemes that probe only the center of the trap, will 
be crucial for the success of the method. 

Further challenges appear when going further towards 
the BCS limit. The obscuring effect of the pseudogap 
and the useful (for this method) effect of the lifetime 
broadening decrease hand in hand, thus in principle there 
is no essential difference to the scenario presented here 
|l6| . However, in the extreme BCS regime only atoms 



close to the Fermi surface are paired and the atoms far 
from it will be transferred resonantly without the shift of 
the peak by the gap energy. In this case, Pauli blocking 
the lowest energy states for atoms |e) by having them 
initially populated would strengthen the signature. 

In summary, we have studied the possibility of observ- 
ing resonance superfluidity by a probing field transferring 
atoms out of the superfluid. The emergence of the pseu- 
dogap obscures the phase transition signature because 
it makes the transfer resonance peak asymmetric even 
above Tc. However, the incoherent nature of the pseu- 
dogap broadens the quasiparticle peaks in the spectral 
densities, allowing one to distinguish between the pseu- 
dogap and superfluid regimes. 
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